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Abstract — In this work we apply the Newton type of local and 
nonlocal internal body fluxes to create the constitutive law for 
the internal body fluxes. We obtain the integral differential 
equations of the stated heat conduction problem and consider 
the steady and unsteady linear and nonlinear one dimensional 
problems. The analytical solutions of the problems are obtained 
using theoretical analysis, differentiation of the integral 
differential equation and the method of separation of variables. 


III. Statement of heat conduction problem 
Consider the classical heat conduction problem with the 
equation according to [1] 


U* 2 


*1 + + l"«o +9i - c o^^7 (!) 


Index Terms —Integral differential equation, local and 
nonlocal internal body fluxes. 


I. INTRODUCTION 

The statement of the problems in heat conduction usually 
includes the surface heat fluxes inside the body and they do 
not consider the constitutive law for internal body flux 
[l]-[4] .The body fluxes are considered as the internal or 
external heat sources. Then the linearized theory must accept 
the nonphysical singularities in temperature field. The 
introduction of the internal body fluxes allows improving of 
the heat problems at least in the sense of excluding the 
nonphysical point singularities. 


II. Internal surface and body fluxes 

Consider a body and let us take some control volume, 
which includes a fixed number of particles. The control 
volume is surrounded by a control surface. The particles 
which are inside the control surface are called internal 
particles and they belong to the control volume. The particles 
which are outside the control volume are the external particles 
and they do not belong to the control volume. All other 
particles belong to the boundary particles of the control 
volume. There are interactions between particles for example, 
according to Newton’s law of cooling; the resultant of 
interactions applied to all internal particles of the control 
volume from the external particles is the internal body flux. 
The interactions applied to the boundary particles of the 
control volume from the external particles are the surface 
fluxes. The Fourier law could be accepted for the internal 
surface heat fluxes and the Newton’s cooling law is taken to 
describe the nonlocal body heat fluxes. 


where x,y,z are the Cartesian orthogonal coordinates, t is 
time, u(x,y,z,t) is the temperature, p(x, y, z) is the 
mass-density of the body per unit volume, Cq is the specific 
heat, k is the coefficient of thermal conduction, q 0 is a rate 
of internal body heat flux per unit volume, q\ is a rate of 
internal heat generation per unit volume produced in the body. 

The introduced in (1) term q 0 could be taken using the 
Newton’s law of cooling in the form of a sum of local and 
nonlocal fluxes 

q^C+gZ 10 ', ( 2 ) 

where 

q'° c = ccpi^x , y , z'), ( 3 ) 

(4) 

= J/J «2 [u{x,y,z,t)-u^.q.C.t^dgdqdC 

where vis the volume of the body, a^,a are constants in 

Newton’s law of cooling, which can depend in general on the 
coordinates of the body. Let us take them constant, p >q is a 
constant The much more general internal body flux can be 
taken into account using the similar approach as in [5]. 

The correspondent initial and boundary conditions should 
be added to create a well posed initial boundary value 
problem. 

It was shown in [2] that 2D and 3D heat problems have 
nonphysical solution in case of a given point boundary 
condition. This forces to look for a new model which could 
include the physical solution into consideration. 
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IV. Methodology and results 
In this paper we consider two simple problems in the 
suggested heat conduction model. The first one is a one 
dimensional steady state problem and the second one is about 
a time dependent solution. 

Note, that the integral term in differential equations causes 
difficulties in analysis. Therefore, we have to modify and 
generalize classical methods. As the result the analytical 
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solutions of the problems are obtained by separation of 
variables and differentiation of the integral differential 
equation. 


V. Problem 1. Steady state case 
Consider the particular case of the problem stated in section 3. 
The equation is 

-a 1 u-^\{uix')-u{,s'))ds = 0 ( ) 

where 0 < x < l . 

The equation (5) can written in the form 

. d 2 u a f . . , _ (6) 

k—--au + —^\u(s)ds = 0, 

dx 2 l p { 

Where 

or, 

a=a ^jk- 

(7) 

The boundary conditions are taken as follows 

u( 0) = w 0 , u(l ) = U\ (8) 

Differentiating (6) one gets 


; +C e +c 


The system of linear algebraic equations (13), (15), (16) 
has the following solution 


A = m\ e ^ - e ^ | + 2| 2 — e ^ 


=(u 1 +u c 


< 19 > 

A 2 - 1 jfu -u a -mu^ + m(u 1 -u Q ), 

(20) f J& \ | 

A 2 = e -1 \lu 1 -u 0 +mu 0 \-m(u 1 -u 0 ). 


The case / — > oo can be obtained from the solution (10) 
and the conditions (8) if C 1 =0,C 3 =0. Then 


ku-cu' = 0 (9) 

The general solution of (9) is 



where Cj, C 2 , C 3 are constants of integration. 

These constants can be obtained satisfying (6) and the 
initial conditions (8). 

Substituting (10) into (6) one gets 



If the length is infinite then (11) will be an identity only 
if 


C 1= 0,C 3 =0. 


( 12 ) 

Consider the case of finite length. If the length is finite 
then the integral in (11) is zero and integrating we obtain 


-mC 1 +C 2 




(13) 


where 



Substituting (10) into the boundary conditions (8) we get 

C 1 +C 2 +C 3 = 0, (15) 


u(x) = C 2 e ^ ( 22 ) 

Using the conditions (8) we obtain the representation of 
the solution (10) of the problem (6)-(8) in the form 

u(x) = u 0 e ^ ( 23 ) 

and 

u % =0. ( 24 ) 

Remark. If Z = oo and a = 0 then we have the classical 
case of the problem and the continuous solution in this 
case does not exist. 


VI. Problem 2. The time dependent problem 
Consider the equation 



subject to the boundary conditions 

u{ 0, t) = u(l,t) = 0 (26) 

and to the initial condition 

u(x,0) = f(x), (27) 
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where fix) is the given distribution of the temperature 
at the time t = 0. 

To obtain the solution of the problem (25)-(27) we apply 
the method of separation of variables which gives us the 
representation of the solution in the form 


u(x,t) = X(x)-T(t) 

Substituting (28) into (25) we get 

cXT' = kX"T - a XT + ^ f XTdx 

1 P 0 

and dividing (29) by XT yields 


(28) 


r 

k X" 

_a_ a 2 (• 

T ' 

~c 0 X 

c* cfx\ 


where A is a constant. 

The boundary conditions (26) imply that 

X(0) = 0, X(Z) = 0 (31) 

To determine A consider (30). 

kX"-aX+°^[xdx = c n AX (32) 

l p { 

Differentiating (32) we obtain 

kX'"-{a + c 0 X)X'= 0 (33) 

Now we consider (33) subject to the conditions (31). 
The following cases are possible: 
a) Let a + c 0 A < 0. 


Denote k 


CC + CqA 

:,= V ~ 


then the solution of (33) is 
(34) 


X = Cj + C 2 cos k]X + C 3 sin k\X 
Substituting the solution (34) into (32) we get 

a f*r -| (35) 

-(o^ +c a X)C 1 + yj [C 2 cosA 1 Jr+C 3 sinA 1 JrJc&r=0 

or 

-mC x + C 2 sin/q/ - C 2 (cosk x l -1) = 0, ( 3( >) 


(a. +cX)l p k 
m =— 1 - 0 - L . 


The conditions (31) imply that 

c x +c 2 = 0 , 


C cosk I+C sink I+C =0. 


Consider the determinant of the system of linear 
homogeneous equations (36), (38), (39). 


-m sin kj -cos^Z + l 
11 0 
1 cos k 1 sin k 1 


The equation (40) yields 


-sin/r./f /?7 + 2tan— 
1 2 


and the roots of (41) can be found analytically 
sinAj/ = 0, k x l = nn x ,n x =0,1,2,... 
n 2 kn f + al 2 

X =- L -- ,n x = 0,1,2,... 

n i ^ /2 1 


and numerically 

k 1 m 

tan—— =-, 

2 2 


kk. + a 

X =-1-,/?= 0,1,2,... 


(44) 

Then we obtain C 2 = —C\ and C 3 is an arbitrary 
constant. 

sin £7 


c=r- 


3 cos k.J-1 


Then the solution of the problem (31), (32) consists of 
two sets of functions 


X n {x) = C n sin^y^,^ =0,1,2,... 


—(l-cosTV A') + sin7V x\, 

\ 


7 


b)Let a + c 0 A = 0. 

Then (33) becomes X m =0 and 
X (x) = Cj + CjX + C^x^ 
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Substituting (49) into (32) we get 

^cM i+ c 2 , + c 3 ^--o ™ 


Ci+ c 2 L +c [L + ^1 

1 2 2 3 3 « 


= 0 . 


(51) 


The boundary conditions (31) imply 

c t = 0, (52) 

C 1 +C 2 I + C 3 J 2 =0. (53) 

The solution of the linear algebraic system of equations 
(51), (52), (53) is 

c x =0,c 2 =0,c 3 =0 ( 54 ) 

and only zero solution exists in this case 

X(jy) = 0. ( 55 ) 

c)Letor + c 0 A > 0 . 

Denote k 2 = a + c °^ j then the general solution of 
V k 

(33) is 

X = Cj + C 2 e klX + C 3 e~ k2X (56) 

Substituting (56) into (32) we get 

_ («, +c a l)l? +c j e ^ dx+c j e -^ dx=Q < 57 > 

a 7 J n : n 


+^C n |^™^l-coSiV n A-j+sinA^ 

The initial condition (27) should be used to obtain the 
constants^ £ . Then (27) and (65) imply that 


f ( X >ll o C n 1 { Sin -j Lx j + 

+'Y J C n j^y^l-cosA^ xj + sinA^ A'j. ( ' 66 ' ) 


The Fourier coefficients are as follows 

nn.x 


c =^~ 


J/(x)sir 


dx 


f . 2 nnx 

J sin— -j—dx 


N x \cL 

02 i 


{/(y)^(i-cosW w) + sinA 7 n x |dx 
if—(I- cosjV wWsinAf x \ dx 


VII. Conclusion 

The main results of this work are the obtained analytical 
explicit forms (10), (65) of the solutions of the new statements 
of the heat conduction problem. 


W , Cj(g -.. 1)= ," (^> 

a 2 

The boundary conditions (31) are 

^+^+^=0, (59) 

C x + C 2 e kl ‘ + C 3 e~ k?l = 0. (60) 

The dete3rminant of the system of linear algebraic 
equations (58), (59), (60) equals 




> 0 . 
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The solution of the linear algebraic system of equations 
(58), (59), (60) is 

c x =0,c 2 = O ,c 3 = O ( 61 ) 

and only zero solution exists in this case 

X<X) = 0. (62) 

Now find the correspondent functions T n (t) using (30) 
and 

A ,A («) 

V A 

according to (42), (44). 

We have T' n - A n T n = 0 and 



T n (t) = A n e ; " ( ■ (64) 

Then the following formal solution of the problem (25) 
-(26) is obtained 

„ . ^ ^ f nn.x'] xt 

K*,o=X c J sin j f 1 +■ 
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